Classification of Critical Phenomena having a Parameter-Dependent Renormalization 
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A generic classification of critical behavior is provided in systems for which the renormalization 
group equations are parameter dependent. It describes phase transitions in a wide class of networks 
with a hierarchical structure but also applies, for instance, to conformal field theories. Although 
these transitions generally do not exhibit universality, three distinct regimes of characteristic crit- 
ical behavior can be discerned that combine an unusual mixture of second-order, infinite-order, 
and discontinuous transitions. In the spirit of Landau, the problem can be reduced to the local 
analysis of a cubic recursion equation, here, for the renormalization group flow of some general- 
ized coupling. Among other insights, this theory explains the often noted prevalence of so-called 
inverted Berezinskii-Kosterlitz-Thouless transitions in complex networks. As a demonstration, a 
one-parameter family of Ising models on hierarchical networks is considered. 



The renormalization group (RG) [2, 0] is by now a 
method found in any "classical" statistical physics text 
book [3-5]. It has allowed to categorize broad classes 
of equilibrium systems into an enumerable set of uni- 
versality classes, each characterized by discrete features, 
such as their dimension and the symmetries adhered to 
by their Hamiltonians. Such universality is made possi- 
ble through the property of scaling that is an inherent 
feature near phase transitions [(| exhibited by these sys- 
tems in regions of the space spanned by their physical 
couplings. Scaling entails that system-specific details on 
the microscopic level become irrelevant, as the behavior 
over many orders in the range of the interactions become 
self-similar. In this framework, analogous behavior in a 
surprisingly wide set of phenomena, such as the conden- 
sation of fluids, spontaneous magnetization of materials, 
or the generation of particle masses in the early universe, 
can be described with a few effective theories - a major 
intellectual accomplishment of modern physics 

Unlike the Euclidean arrangement of atoms in a lattice, 
agents in biological or social systems may exhibit com- 
plex networks of mutual interactions ■ The study of 
critical phenomena is inseparable from the understand- 
ing of the geometry of the network [lol ]. It has been 
realized that many of the networks that are engineered 
by some natural or human activity themselves exhibit 
emergent complex properties, for instance, as found in 
the scale- free degree distribution of the internet. While 
these networks, and dynamical systems on them, may 
behave critically, those phenomena were soon found to 
be non- universal i-e., they are intimately tied to in- 
tricate details of the specific system. In this sense, it 
would seem unlikely that a sweeping classification could 
be devised. In this Letter, we will categorize equilib- 
rium phenomena observed on a large set of networks 
having hierarchical structure. Our discussion explains, 
for example, the robustness of infinite-order transitions 
in a range of quite distinct network models summarized 
in Ref. [9J, or in conformal field theory, where it sig- 
nals the loss of conformality [12]. Our study shows that 
criticality in these models is generally non-universal but 



falls into three generic regimes, here characterized by 
their divergence in the correlation length [l^. One of 
these regimes is an infinite-order transition reminiscent 
of that described by Berezinskii and Kosterlitz and Thou- 
less (BKT) 0| . It is flanked on one side by a transition 
with a weaker, algebraic divergence, similar to a second 
order transition (albeit non- universal) , and on the other 
by a regime with an even stronger essential singularity. 

Our classification immediately applies to hierarchical 
networks on which RG is exact and transitions can be 
studied in great detail. There, these regimes are char- 
acterized by the relative strength of long-range, small- 
world bonds 0]- Metric version of such networks, like 
the Migdal-Kadanoff RG, provide classical text-book ex- 
amples of RG and universality 0, [HI, 15 1. But in the 
advent of complex networks, many hierarchical designs 
with non- metric (i.e., small- world or scale- free) proper- 
ties have been devised and studied flB - El! ]. 

The central tenant of real-space RG consists of a pro- 
cedure whereby the partition function of the original sys- 
tem is mapped recursively onto itself after tracing out a 
fraction 1 — 1/6 of the dynamic variables, in some form 
of "blocking". Prior couplings 4 between those variables 
combine non-trivially to produce new, effective couplings 
K n +i between the remaining variables after the n-th RG- 
step while leaving the Hamiltonian form-invariant. This 
mapping constitutes the RG-flow 



-i =U («„) . 
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where 1Z indicates a (typically non-linear) set of recur- 



sions. In the thermodynamic limit, 



\og b N -> oo, 



phase transitions are characterized purely by a local anal- 
ysis for K n ~ K n +i ~ near the fixed points (FP) of 
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independent of kq. Here, ko represents the "raw" cou- 
plings of the original system that carry the dependence 
on the system's control parameter /i £ [0,1], referring to 
a temperature, Ko 0-0 = fi = e~P J in units of J = 1, 
or a bond-percolation probability, po = p = 1 — say. 
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Since kq expresses microscopic details of a potentially 
large family of conceivable systems adhering to Eq. (fl}, 
the insensitivity of k^ on k is an expression of univer- 
sality: only certain symmetry properties of the original 
systems remain preserved by 1Z. Therefore, the linearized 
expansion, K n ~ Koo + e*n with small e„ for large n, near 
the FP provides a full accounting of the macroscopically 
observable properties of any such system via the eigen- 
value problem obtained from 
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The eigenvalues A of the Jacobian ^5 (koo) and their 
eigenvectors u\ provide the scaling exponents and scaling 
fields observed in the phase transition. 

A much richer phenomenology arises when the RG-fiow 
1Z itself becomes dependent on the control-parameter. In 
that case, Eq. (fTJ) generalizes to k„+i = 1Z (k„; fi) with 



K oo — 7l> (^oo, A*) 



(4) 



i.e., the FP k^ becomes a nontrivial function of fi. The 
consequences of such behavior (for a single coupling) are 
depicted in Figs. [1] First, consider the case of constant 
FP shown in Fig. HJa) . Drawing constant FP as function 
of the control-parameter may seem redundant, however, 
it allows to illustrate the connection to the raw couplings 
Kq (fi) (green-dashed line), which must rise smoothly with 
fi from k (0) = to k (1) = 1 to be physical, but are 
otherwise arbitrary (here, chosen as a straight line). Be- 
low (above) the point where kq (fi) intersects the unstable 
FP, the RG-fiow evolves toward the stable FP at k^ = 
(koc = 1)- The arbitrariness in kq (fi) reflects the univer- 
sality in the family of systems obeying the same FP: no 
matter at which value of fi c the specific system kq (fi) 
intersects the unstable FP with < k^ < 1, its (it- 
independence guarantees identical critical behavior. 

We claim that the remaining three panels in Fig.[T]cap- 
ture all generic features that can arise for /i-dependent 
FP. Variable FP Koq (fi) can collide, either linearly or in a 
square-root branch point (BP); any other behavior (such 
as a higher-order BP) is exceptional. We can devise a 
simple theory[26] that reproduces these generic features, 
thereby demonstrating the generality of this classifica- 
tion, not only accounting for hierarchical networks but for 
any physical system described b y p arameter-dependent 
renormalization group equations 12, HI, TB- 21|. 

It is sufficient to consider the RG recursion for a single 
coupling n n with some control parameter fi. We argue 
that FP in a real, /i-dependent RG-fiow 1Z in Eq. ([J} 
will exhibit BP at some point(s) (fi B ,K B ). Near fis we 
express generically 1Z (k; fi) ~ a (fi) k + b (fi) k 2 + c (fi) k 3 
because the need for a strong-coupling solution = 
prevents a constant term and requires at least a cubic 
form to achieve BP. With generically analytic coefficients 
at fig, we expect a (fi) ~ ao + a% (fi — fis), b (fi) ~ bo, 
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Figure 1: Generic plots of the fixed points (FP) 
function of fi (a) for conventional FP, (b, c) for FP with a 
physical branch point (BP), and (d) for FP with BP outside 
the physical domain (0 < fi, n < 1). In (a), /i- independence of 
FP ensures universal critical behavior wherever the raw cou- 
plings kq (fi) (green-dashed line) intersects the unstable FP 
(red line); a blue dot and dashed line mark the critical point 
fi c , and blue arrows indicate the RG-fiow from Ko (fi) toward 
the nearest stable FP (black lines). In (b), Ko (fi) still inter- 
sects at an unstable FP below BP, such that the RG-fiow does 
not pass BP, leading to quasi-conventional behavior but with 
/i c -dependent exponents. In (c), ko (fi) is located above BP so 
that the RG-fiow must pass near BP at fi c , leading to BKT- 
like behavior. In (d), BP drops below the physical domain 
and only stable FP are accessible, resulting in an exponential 
divergence at the intersection fi c of two FP-branches. 



c(fi) ~ Co for fi — > fiB- Locating BP at (/is,Ks) fixes 
ao = 1 + cqk b and bo = —2coKb- To orient BP correctly 
requires ai/co < 0, and we set a\ — —cqA 2 with A > 
0. Finally, stability of the strong-coupling FP at 
demands cq < 0, and we may set cq = — 1. This yields 
the minimal model 

K n+ i-K„ ~ ^ ~ [-«! + A 2 (fl - fl B )] K n +2K B Kn-K S l . 

An 

(5) 

After extracting = [27j , the remaining FP equation 
indeed produces a BP at (fiB, &b)'- 



k b ± A^Jfi - fi B . 



(G) 



Local expansion near each FP as in Eq. ^ provides the 
eigenvalues A (fi) = 1Z' («oo; m) depicted in Fig. [Ha), 

A = 1 - k% + A 2 (fi - fi B ) , 

X ± = 1 t 2Ak b - Mb - 2A 2 (ji - fi B ) , (7) 

In Fig.[TJ panels (b) and (c) correspond to cases where 
BP at (fiB, Kb) is in the physical domain (0 < /i, k < 1); 
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Figure 2: Plot of the eigenvalues (a) in Eqs. ([7J for the RG- 
flow in Eq. JHJ, here A — 1, kb ~ §, and (b) for the RG-flow 
in Ref. fjj] for HN5 (at y = 0.1). At BP (fj, = u B ), two 
conjugate eigenvalues emerge simultaneously with marginal 
stability, A (pb) = 1, and A - remains unstable (> 1, marked 
red) until the lower FP-branch, drops below k^. In (a) 
for kb < 0, setting kb — > -kb merely swaps A* — ¥ A^ 
(see Eq. ([2])), and now A (fj, c ) = A + (p c ) = 1 intersect at 
Mc > Ms with marginal stability. For HN5 in (b), A* also swap 
positions when kb becomes negative but in a less symmetric 
manner (see Supplement B). 

panel (d) represents kb < 0. Within this domain, the 
lower FP-branch (/i) is always unstable (A - > 1, see 
Fig. [2]) while (m) remains stable. Stable and unsta- 
ble branches merge at BP, where particularly interesting 
phenomena arise. The decisive difference between panels 
(b) and (c) is the location of BP relative to kq (m). 

For the case of panel (b) (ej^ when long-range, hierar- 
chical couplings are weakest [13]), Mb is small and/or Kb 
is closer to unity (or even above) . Then, kq (m) merely in- 
tersects the unstable branch (m) at some critical point 
Me > Ms defined by (/i c ) = (Mc)- The RG-flow 
(vertical blue arrows in Fig. [IJ for < [i < Mc advances 
toward strong coupling, k^, while for \i c < \i < 1 it 
flows toward /c+ (m) [28|. Near Mo the critical dynamics 
of the system is now determined by the local properties 
of the unstable FP (m c ) that has been selected by the 
specific system via its raw coupling kq (m) • As for a con- 
ventional system in Eq. © , local analysis [5] of Eq. ([5]) 
near (m c ) yields the diverging correlation length, 

£~|M-McP (Pc) , M^Mo (8) 

but with a non-universal exponent y t — log 2 A~ (/i c ) = 
l/V(Mc)- For /ic \ Mb ; ^ becomes marginal and the 
exponent diverges as v (/i c ) ~ l/VMe ~ Mb- Yet, for [i c > 
Mb, the RG-flow never passes sufficiently near BP to be 
affected. 

For the case of panel (c) where kq (/i) passes above BP 
(e.g., for somewhat stronger long-range couplings [13]), 
the RG-flow must pass BP which now dominates crit- 
ically, i.e., Mc = Mb, with an infinite-order divergence 
characterizing this regime. Well below (above) mb, the 
RG-flow evolves unperturbed to (to k+ ), the closest 
stable FP. However, just below mb the RG-flow gets ever 
more impeded near BP before it can reach k^. Asymp- 



totically for m /* Mb near K„~KB + e„ with small e„ at 
large but intermediate n, Eq. ((5]) provides 

£n+i - e n = ~ -K B A 2 (mb - m) - K B € 2 n . (9) 
An 

This relation exhibits a boundary layer, i.e., in the limit 
M /*" Mb the solution drastically changes behavior. With 
the methods of Ref. [22] , we rescale e„ — > 7e„ and n — > Sn 
to obtain a balance for S ~ 1/7 ~ 1/Vmb — M- Ac- 
cordingly, the characteristic width of the boundary layer 
scales with n* ~ 1/ %/Mb — Mj which leads to the diver- 
gence in the correlation length characteristic of BKT, 



£(M)~2" ~ e V^, m^Mb=Mc- (10) 

Clearly, the physical origin of this singularity is not re- 
lated to an actual BKT transition, with its formation of 
delicate topological structures In fact, instead of be- 
ing rare, it is one of three generic types of transition often 
found in hierarchical networks [9]. 

The most unconventional behavior is depicted in panel 
(d) of Fig. m where kb < (e.g., when long-range cou- 
plings dominate [13j). No unstable FP can be reached 
for any physical choice of kq (m). The RG-flow always 
advances to the closest stable FP, either at strong cou- 
pling, kP^ for < m < Mc; or at patchy order, k+ (m) 
for /i c < /i < 1. Both lines of FP cross at Mc(> Ms), 
defined by k+ (m c ) = kP^ = 0. This condition im- 
plies that both their eigenvalues are equal and marginal, 
A (Mc) = ^ + (Mc) = 1) as must invert its stability at 
the intersection, making marginal stability inherent to 
any such system. In our model, k+ (m c ) = in Eq. (|6|) 
provides —kb = y//J. c — Ms, hence, Eqs. ([7]) give A° :+ ~ 
I- A 2 |m - Mc| for M -> Mc see Fig.^a). Near k°^+ (m c ) = 
0, the local analysis on Eq. according to Eq. ([3]) gives 
~ A°'+e„ or e„ — e exp (-nA 2 |m - Me|) with a 
cross-over at ?i* ~ l/\fj, — [i c \ that is generic when A 
and A + intersect linearly. Thus, the correlation length 
diverges as 

£(m)~2"* = e ra, n^nc (11) 

Again, the RG flow does not pass BP, since it is located 
below the physical domain. 

We can extend the discussion to include the effect 
of further control parameters, such as an external field 
rjo = r) = e _/3 ' 1 . We find that its generic RG-flow can be 
expressed asymptotically as 

7/ n+1 ~r^, X h ~ 2 - C K + (mc) (12) 

with some Mc-dependent constant C > 0, near the criti- 
cal point Mc and for sufficiently small k^ (mc)- (E.g., in 
Supplement B we find C = 2m c for the Ising model on 
HN5.) This satisfies the physical requirements on its FP, 
iloo = 0, 1; only for h — the RG-flow remains at the 
unstable FP, i.e., r\ n = 1 f. a. n, and for any h > the 
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stable strong-coupling FP at rj^ = is reached eventu- 
ally. The eigenvalue near the unstable FP satisfies A^ < 2 
such that yh = log 2 A/j < 1, as shown in Ref. (23j. There, 
a scaling theory is developed concerning the BKT regime 
based on the exponentially divergent correlation length in 
Eq. (|10p . leading to an order-parameter (magnetization, 
fraction of sites on percolating cluster, etc.) 

t-i / const. (1 - y h )\ 
m ~ £ - exp <^ — — \ , n / fj, c , (13) 

I (A*c-M) J 
when y t — > 0. Of course, for y t > it is [f| 

o 1 - 



m ~ (/x c - /«) 



.'7/ 



(14) 



for small- world systems where AT takes the role of L d [23| . 

Our theory in Eq. ([5]) not only explains the robust- 
ness of the value of Xt = i conjectured in Ref. (23| . 
but also broadens the scope to a total of three generic 
regimes in the divergence of £, as we have explained. 
With the addition of Eq. (fT2")l , we can account for the 
behavior of the order-parameter to. For the first two 
regimes where Kb > 0, it is Ah < 2 in Eq. ([TO so that 
1 — y/j > 0. In the regime with the weakest distortion of 
the FP, see Fig. [ljb), we have shown in Eq. ((HJ that 
yt = l/^(/i c ) > 0, i.e., it is < B (p c ) < oo simi- 
lar to an ordinary 2nd-order transition, except for its 
non-universal /^-dependence. In the BKT-regime, see 
Fig. [He), it is fi c = fj,B and A + (/i c ) = 1 in Eq. (J7J such 
that yt = a nd B — > oo, which leads to Eq. (fl~3|) de- 
scribed in Ref. [23]. Finally, when kb < 0, see Fig. [lid), 
it is n c > hb and k+ (// c ) = 0. And although A + (p c ) = 1 
leads to the divergent correlation length in Eq. ([TO , it is 
A/,. = 2 in Eq. (fTSj) and B = 0, hence, the order-parameter 
to in Eq. (|14jl becomes discontinuous. Such an unusual 
discontinuity, originating in intersecting stable FP, has 
indeed been shown for ordinary percolation on a hier- 
archical network 24| . We plot the magnetization of the 
Ising model on a hierarchical network, which follows 2nd- 
order behavior, Eq. (fT4f . in the first regime in Fig. [3][a), 
it has an infinite-order transition, Eq. (|13p . in the BKT 
regime in Fig.[3][b), and it develops a discontinuity in the 
regime of intersecting stable FP in Fig. [3tc) . 

As a demonstration for our theory, we revisit the Ising 
model on the Hanoi network HN5 previously considered 
for h = in Ref. [13]. There, a one-parameter fam- 
ily of Ising models was conceived via the ratio between 
short-range and small-world coupling strengths that in- 
terpolates between all three regimes; Fig. 11 in Ref. [l3j 
corresponds to Fig. [T][b-d) here. That system is far more 
complex than our model here in that there are two cou- 
plings and three fields (when h > 0, for site-, bond-, and 
three-point magnetizations) to be renormalized. Yet, the 
same three regimes in the divergence of the correlation 
length £ and the magnetization m ensues, as our theory 
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Figure 3: Plot of the magnetization in the Ising model on the 
hierarchical network HN5 from Ref. [13J (a) for y — 0.1, (b) 
for y = 0.4, and (c) for y = 1, which in sequence correspond 
to regimes (b) to (d) in Figs. fH as shown in the supplement. 
The transition is (a) 2nd-order with /3 = 0.205 . . ., (b) BKT- 
like, and (c) discontinuous. Finite-size effects (N ~ 2 fe < oo) 
remain throughout. 



predicts. In Supplement B, a Mathematica script repro- 
duces plots identical to Figs. Q] and [5] for this case; more 
examples will be presented elsewhere. 

Conclusions: Ref. (23| has raised the question about 
the universality of the BKT result in Eq. (fTU|) . Our dis- 
cussion shows that the robustness of this behavior derives 
from the fact that BP is most generically a square-root 
singularity. Only if for some reason the cubic and higher- 
order terms in the RG-flow near BP in Eq. ([SJ were ab- 
sent, a higher order root-singularity could arise. The full 
exponential singularity in Eq. (jlip is even more solid, as 
it does not depend on the nature of BP but merely in the 
fact that two intersecting lines of FP must switch stabil- 
ity, which implies marginal eigenvalues at the point of 
collision. While those eigenvalues scale linearly with fi 
there, the exponential divergence in Eq. (fTTj) is robust. 
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